A study on the anomaly of p over n ratios in Au + Au 
collisions with jet quenching 



Xiaofang Chen 1,2 , Hanzhong Zhang 1,2 ' 3 , Ben- Wei Zhang 1,2 and 
Enke Wang 1 ' 2 

institute of Particle Physics, Huazhong Normal University, Wuhan 430079, China 
2 Key Laboratory of Quark & Lepton Physics (Huzhong Normal University), Ministry 
of Education, China 

3 Department of Physics, Shandong University, Jinan 250100, China 

Abstract. The ratios of p/n at large transverse momentum in central Au + Au 
collisions at RHIC are studied in the framework of jet quenching based on a next-to- 
leading order pQCD parton model. It is shown that theoretical calculations with a 
gluon energy loss larger than the quark energy loss will naturally lead to a smaller p/n 
ratios at large transverse momentum in Au + Au collisions than those in p+p collisions 
at the same energy. Scenarios with equal energy losses for gluons and quarks and a 
strong jet conversion are both explored and it is demonstrated in both scenarios p/n 
ratios at high pt in central Au + Au collisions are enhanced and the calculated ratios 
of protons over pions approach to the experimental measurements. However, p/p in 
the latter scenario is found to fit data better than that in the former scenario. 
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1. Introduction 

Relativistic heavy-ion collisions provide the possibility of creating in the laboratory a 
plasma of deconfined quarks and gluons, or QGP. So far a large amount of interesting 
data have been accumulated in the experiments at the Relativistic Heavy Ion Collider 
(RHIC) [H El [31 H], which strongly suggest a new kind of matter may have been 
formed at RHIC. One of the most striking evidences of the formation of QGP is the 
strong suppression of hadron spectra with large transverse momentum (or high p?) in 
Au + Au collisions as compared to p + p collisions [51 [6]. This kind of suppression 
agrees with the prediction of jet quenching theory [7J El El [Hll E], which has provided 
compelling explanations on many novel phenomena at high p? observed at RHIC, such 
as the suppression of 7r° production at high p T , and the disappearance of back-to-back 
correlations of high pt hadrons [T2] . 

In the picture of jet quenching or parton energy loss, an energetic parton will loss 
energy by multiple scattering in the medium when it propagates through the hot QGP 
medium. The majority of current approaches to the medium-induced parton energy loss 
may be divided into four major schemes which included higher twist (HT) [TB I [T4l [TBI [16] . 
path integral approach to the opacity expansion (BDMPS-Z/ASW) pH [TBI UHl EQ] , finite 
temperature field theory approach(AMY) [211 1221 [23] and reaction operator approach 
to the opacity expansion (GLV) pU ES, EHl ETj. It has been shown that the energy 
loss of light quarks and gluons is proportional to the gluon density and has a quadratic 
dependence on the total distance traversed by the propagating parton due to the non- 
Abelian LPM interference effect in multiple scatterings [TBI EH EHl [291 SO] • Moreover it 
is expected that the gluon energy loss in nuclei is 9/4 times of the quark energy loss due 
to the different color factors for quark-gluon vertex and gluon-gluon vertex [3 EHl 30J. 
Because high p? proton (anti-proton) production in p + p collisions is dominated by 
gluon fragmentation whereas high pt pion at p + p collisions is dominated by quark 
fragmentation, it is expected that larger gluon energy loss than quark energy loss in 
hot /dense medium will lead to a smaller p/ir ratio at high px in Au + Au collisions than 
in p + p collisions at the same energy [31] . 

However, STAR Collaboration recently published very interesting measurements 
which indicate that p/n + and p/tt~ ratios at high pt in central Au+Au collisions [32] 
approach those in p + p and d + Au collisions [33], and the nuclear modification factor 
of protons is similar to that of pions [32]. To solve the discrepancy between the 
experimental data and the theoretical calculations, a jet conversion mechanism [34j has 
been proposed by taking into account the possibility of a quark jet in QGP converted into 
a gluon jet and vice versa, which has previously been considered theoretically in Refs. 
[35] for eA deeply inelastic scattering within higher twist expansion approach. It is found 
that conversions between quark and gluon jets indeed lead to an increase in the final 
number of gluon jets in central heavy ion collisions than the case without conversions, 
but to explain experimental data observed by STAR a conversion enhancement factor 
K = 4 ~ 6 would be needed [34] . 
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In this paper, we will study the p/it ratio in Au + Au collisions with yfs = 200 GeV 
by using a next-to-leading order (NLO) pQCD inspired parton model, which has been 
very successful in describing the high px region physics in p + p collisions [36] . Later by 
including the nuclear effects it is extended to consider the single hadron and di-hadron 
productions in Au + Au collisions and has given a very good description on the high pt 
pion and charged hadron productions at RHIC [37J. An interesting observation in the 
model shows that the ratio of gluon/quark jets in NLO is larger than in LO calculations 
[37]. In the same framework, we will investigate the high px pion and proton (anit- 
proton) spectra, especially the p/ir ratio at yfs = 200GeV. We find that with gluon 
energy loss larger than quark energy loss in QGP, the NLO pQCD parton model shows 
the ratios p/n in Au + Au will be smaller than those in p+p collisions, which agrees with 
the previous pioneer study in leading-order parton model [31]. To explain the observed 
p/tt ratios at STAR Collaboration phenomenologically, two different scenarios in jet 
quenching theory are considered: one is the case that gluon energy loss is equal to quark 
energy loss in nuclear medium; another with a strong jet conversion mechanism where 
a appreciable quark jets are converted into gluon jets. Both scenarios will effectively 
increase the total number of final gluon jets passing through the hot QGP medium, 
and thus give a higher p/n ratios, which are found to reproduce the experimental data 
fairly well. However, p/p in the latter scenario is found to fit data better than that 
in the former scenario. In addition, RP^ = R\a * s found in the former scenario while 
R AA < R\ A m the latter scenario. 

The paper is organized as follows. In Sec. [2, we describe explicitly the NLO pQCD 
parton model in p + p and Au + Au collisions. The effect due to equal energy losses 
of the quark and gluon is taken into account in Sec. [3j In Sec . HI the mechanism of 
flavor conversion between quark and gluon jet is considered. A reasonable conversion 
rate with which the results are in good agreement with experimental data is fixed on. 
Finally, we conclude in Sec. [5] with a summary of present work. 



2. High pt hadron production given by a NLO pQCD parton model 

PQCD parton model has been shown to work well for large pt particle production in 
high energy nucleon-nucleon collisions [38J. The factorization theorem demonstrates that 
the inclusive particle production cross section in p + p collisions can be expressed as a 
convolution of parton distribution functions inside the hadron, elementary parton-parton 
scattering cross sections and parton fragmentation functions, 

da h r 

= Yl J dx a^bdz c{dfi) f a i v {x ai [i 2 )f h i v {x h ,[i 2 ) 

abcdie) 

x — — — y{ab -> cd(e))D h/c{d e) (z c u e - ) , /i 2 ), (1) 
2x a XbS 

where ^>(ab — > cd(e)) is related to the squared matrix elements for the various 2^2 
or 2 — > 3 subprocesses in parton-parton hard scattering. dT is the differential two- 
body or three-body phase space element. The detailed expressions and introductions 
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Figure 1. (Color online) The hadron spectra for (tt + + 7r~)/2 and (p + p)/2 and the 
ratios of (p + p)/(n + + 7r~) in p + p collisions at yfs = 200 GeV. Different values of 
the scales for factorization, renormalization and fragmentation are used in numerical 
calculations. The data are from Ref[32j. 

for ^(ab — > cd(e)) and dT can be found in references [36, 39J. f a / p (x a ,fi 2 ) is the normal 
parton distribution function (PDF) for which we will use the CTEQ6M parametrization 
[40J, and x a is the fraction of the hadron's momentum carried by the parton. D\, c {z c , /i 2 ) 
is the fragmentation function (FF) of parton c into hadron h, and z c is the momentum 
fraction of a parton jet carried by a produced hadron. We will use mainly the updated 
AKK parametrization [JT] for jet fragmentation, which has recently been improved from 
both new theoretical input and RHIC data, especially for 7r ± , K , p/p, K° s and A/ A 
particles. 

The calculations discussed in this paper are carried out within a NLO Monte Carlo 
based program which uses a variant of the phase space slicing technique (36J, H2] . This 
allow the same kinematic cuts in the extraction of the data to be imposed on the 
theoretical predictions. For the lowest-order contributions of 2 — > 2 tree level (a 2 order), 
the NLO correction contributions {o? s order) include 1-loop correction contributions 
to 2 —>■ 2 tree level and 2^3 tree level contributions, which utilize two-cutoff 
parameters |36j, 5 S and 5 C , to isolate the soft and collinear divergences in the squared 
matrix elements. The results in a set of two-body and three-body weights depend on the 
cut-offs, but this dependence cancels when the weights are combined in the calculation 
of physical observations. 

For the calculations of the cross section at fixed values of the transverse momentum 
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Figure 2. (Color online) The spectra of (ir + +ir~)/2 and (p+p)/2 in p+p and Au+Au 
collisions, and the suppression factors in central Au + Au collisions at y/s = 200GeV. 
The data are from Ref. 32 . 



Pt of the produced hadron, we will choose the factorization scale, the renormalization 
scale and the fragmentation scale \i to be all proportional to px- As pointed out in 
Ref. [37] . the calculated single inclusive n° spectra in the NLO pQCD parton model in 
p + p collisions agree well with the experimental data at the RHIC energy with the scale 
in the range \x = 0.9 ~ l.bpx- Shown in Fig. we give the spectra for charged pion 
and (anti-)proton in p + p collisions and their ratios {p + p)/ (vr + + 7r~) as functions of 
Pt- The grey bands are with scales /i = 0.5 — 2.0pT- Not only the NLO spectra but also 
the ratios in p + p collisions are not very sensitive to the scale. So we will choose the 
same scale \x = 1.2pr in the following calculations for both p + p and Au + Au collisions. 

Assuming the N + N cross section can be extrapolated to A + A collisions, the 
inclusive invariant differential cross section for producing a hadron h in a collision of 
nucleuses of types A and B can be written as, 

da h AA ^ 1 



^2 — J d(pdb 2 dr 2 dx a dxi,dz c ^,e) 



abcd(e) 

x t A (r)t A (\r-b\) 

x fa/A(x a , /i 2 , r)f b/A (x b , /i 2 , |r - b|) 

x ^f(ab -> cd(e)) 

2x a Xf)S 

X Ai/c(<i,e)(2c(d,e)) /i 2 , AE c ^ e )) , (2) 

where <ft is the azimuthal angle with respect to the reaction plane and describes the 
direction a final parton jet is emitted in the overlap region between two colliding nuclei. 
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The function t A (r) = -^jpy 1 — r 2 /R 2 is the nuclear thickness function in a hard-sphere 
geometry model, normalized to / d 2 rt A (r) = A which is the nucleon number or mass 
number in a nuclei with radius R = 1A2A 1 ^ 3 fm. f a /A(x a , p 2 ,r) is assumed to be 
factorizable into the parton distribution in a free nucleon f b / N (x,fi 2 ) and the nuclear 
shadowing factor S a / A (x,r) given by HIJING parameterization |43j . 

f a/A (x,fi 2 ,r) = S a/A (x, r)[— / a/p (x,/i 2 ) 

+ (l-§)/a/n(x,/X 2 )], (3) 

where we have explicitly taken into account the isospin of the nucleus by considering the 
parton distributions of a neutron which are obtained from that of a proton by isospin 
symmetry. 

As assumed in Ref. [28l [37], the effect of final-state interaction between produced 
parton and the bulk medium can be described by the effective medium-modified FF's, 



D h/c (z c ,AE c ,n 2 ) = (l-e~<*> 



+ e-<x>L>° /c (^ 2 ), (4) 

where z' c = Pt/(pt c — AE C ), z' g = (L/ X)pr/ AE C are the rescaled momentum fractions, 
AE C is the average radiative parton energy loss and (L/X) is the number of scatterings 
along the parton propagating path, 

POO rjq- 

(L/X)= / -— p 9 (r,b,r + nr). (5) 

Neglecting transverse expansion, the gluon density distribution in a 1-d expanding 
medium in A + A collisions at impact-parameter b is assumed to be proportional to the 
transverse profile of participant nucleons, 

P 9 (r, b, r) = m^[t A { T ) + t A (\b - r|)]. (6) 

According to recent theoretical studies [HI [TU 03] the total parton energy loss in 
a finite and expanding medium can be approximated as a path integral, 

A f] poo t -j- 

AE^(—) ld dr %(r,b,r + nr), (7) 

aL J T r po 

for a parton produced at a transverse position r and traveling along the direction n. 
(dE/dL)id is the average parton energy loss per unit length in a 1-d expanding medium 
with an initial uniform gluon density po at a formation time To for the medium gluons. 
The energy dependence of the energy loss is parameterized as 
dE 

(jfjid = e (£M> - 1.6) L2 /(7.5 + £/p ), (8) 

from the numerical results in Ref. |Hj in which thermal gluon absorption is also taken 
into account. The parameter e should be proportional to the initial gluon density p . 
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Figure 3. (Color online) The fractions of contributions to the charged it and (anti- 
proton from quark and gluon jets with AKK fragmentation functions and DSS 
fragmentation functions in p + p collisions at \fs = 200 GeV, respectively. 



A simultaneous fit to the single and dihadron data constrains the energy loss parameter 
within a narrow range: eo = 1.6 — 2.1 GeV/fm in Ref.|37J. For the following calculations 
we will choose the same parameters as in the reference, p,o = 1-5 GeV, eoAo = 0.5 GeV 
and To = 0.2 fm in Au + Au collisions at ^ = 200 GeV. 

Shown in Fig. (121) are the charged n and (anti-)proton spectra and the suppression 
factors in central Au + Au collisions at y/s = 200 GeV. The suppression factor (or 
nuclear modification factor) which mainly manifests jet quenching effect in our current 
studies is defined as [33], 

R AA = da AA/ d y d2 PT , g x 
(Nunary) daftp/ dyd?p T ' 

where Nu nary is the average number of geometrical binary collisions at a given range of 
impact parameters, (Nb inary ) = J d 2 bd 2 rtA(r)tA(\h — r|). 

Non-Abelian feature of QCD show that gluon energy loss is 9/4 times of quark 
energy loss in the dense matter due to the difference of color factors for gluon-gluon 
vertex and quark-gluon vertex [7J [29], [30] . With the parton energy loss, AE g = 9/4AE q , 
and the energy loss parameter chosen as eo = 1.68 GeV/fm, the calculated large 
transverse momentum (tt + + 7r~)/2 spectra in central Au + Au collisions fit data well 
and the large transverse momentum (p + p)/2 spectra fit data bad in Fig. Also 
shown in the figure are the suppression factors for (n + + 7r~)/2 and (p + p)/2 spectra 
in central Au + Au collisions, respectively. The former is found larger than the latter. 

Shown in Fig. ([3]) are the fractions of the contributions to the charged pion and (anti- 
proton from quark and gluon jets with different fragmentation functions, respectively. 
In the two plots the red solid curves and the red dashing curves are with AKK 
fragmentation functions [31] while the blue dotted curves and the blue dashed-dotted 
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Figure 4. (Color online) The ratios of (p + p) / (ir + + tt~) in p + p and 0-12% Au + Au 
collisions at y/s = 200 GeV. The energy loss parameter is chosen as eo = 1-68 GeV/fm. 
The gluon energy loss is considered as about 2 times of the quark energy loss. 

curves are with DSS fragmentation functions [46J. With the two sets of fragmentation 
functions, our NLO numerical results show that high (anti-)proton is dominated by 
gluon fragmentation whereas high px charged pion is dominated by quark fragmentation. 
It is consistent with a previous study based on Monte Carlo simulations [31 J. 

Because the gluon energy loss is considered as 9/4 times of the quark energy loss, 
the contribution from gluon jets is suppressed more greatly than the contribution from 
quark jets in central Au + Au collisions. The non-Abelian feature of parton energy loss 
causes to a smaller suppression factor Raa for large transverse momentum (p + p)/2 
than that for (tt + + ir~)/2, as shown in Fig. (jSJ). So a smaller p/ir ratio at high p? in 
central Au + Au collisions than that in p + p collisions is obtained, as shown in Fig. 
(TJJ. These conclusions had been previously pointed out in Ref. [31] within LO pQCD 
parton model. 

Since the number ratio of gluon/quark jets created in hard scattering processes in 
NLO calculations is larger than in LO calculations [37], the contribution fraction from 
gluon jets will be enhanced in NLO compared with that in LO. Therefore, the NLO p/it 
ratio is larger than the LO ratio in both p + p and central Au + Au collisions. However, 
even in NLO calculations the p/n ratio in central Au + Au collisions is still smaller than 
that in p + p collisions, as shown in Fig.fjlj), while recent STAR data indicate that p/ii + 
and p/tt~ ratios at high px in central Au + Au collisions [32] approach those in p + p and 
d + Au collisions [33] . To solve this discrepancy between the theoretical calculations and 
experiment measurements, we will explore two scenarios: one is that gluon jets lose the 
same amount of energy as quark jets; another is jet conversion between gluon and quark 
jets. In the following we demonstrate that how p/n and p/p ratios will be modified in 
these two scenarios, and their comparisons with the experimental data. 
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Figure 5. (Color online) The charged 7r meson spectra (upper) and the suppression 
factors (lower) with 2 sets of parton energy loss (AE g = 9/4AE q and AE g = AE q ) in 
central Au + Au collisions at y^s = 200 GeV. The data points are from [331147] , 




Figure 6. (Color online) The separated p and p spectra and the suppression factors 
in central Au + Au collisions at s/s = 200 GeV. The data are from [33l [47] . 



3. Scenario I: gluon jet energy loss is equal to quark jet energy loss 

To explain phenomeno logically the observed p/n ratios given by STAR Collaboration, 
we first consider a case in jet quenching theory: equal energy loss for both gluon and 
quark jets. Though currently there is no such a theory of jet quenching to justify that 
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gluons will lose the same amount of energy as quarks, this toy model could shed light 
on the flavor dependence of parton energy loss in nuclear. We note such a scenario of 
parton energy loss had been explored in a study of the sensitivity of hadron spectra 
suppression to the non-Abelian parton energy loss |48j . 

The hadron yield will be less suppressed with parton energy loss AE g = AE q than 
that with parton energy loss AE g = ^AE q in central Au + Au collisions, so one should 
enlarge the energy loss parameter e to fit data. Shown in Fig. ([5]) are the numerical 
results with 2 sets of parton energy loss (AE g = \AE q and AE g = AE q ) for charged 
pion meson spectra and the corresponding suppression factors compared with data. 
From the suppression factors in Fig. (J5j), it is clear that eo = 1.98 GeV/fm works well 
for fitting data. The parameter value will be used to give (anti-)proton spectra with 
parton energy loss AE g = AE q in central Au + Au collisions. 

AKK parametrization of fragmentation functions just give parton fragmentation 
functions to (n + + 7r~)/2 or (p + p)/2. To get ratios p/n + , p/^~ and p/p we need 
the parton fragmentation function to tt + , tt~, p and p respectively. Theoretical and 
experimental studies [321 [33J [3U H9] show that the fragmentation of quark and anti-quark 
jets produce mainly protons and anti-protons, respectively, and gluon jets contribute to 
protons and anti-protons equally. As in Ref . [33l 131] , we assume that no anti-protons 
are produced from a quark jet and no protons are produced from an anti-quark jet. We 
further assume that high pt charged ir + and ir~ mesons are equally produced from hard 
parton jet fragmentations, consistent with experiment findings [32| [33]. 

Based on the NLO pQCD parton model, we calculate the separated proton and 
anti-proton spectra with the energy loss parameter eo = 1.98 GeV/fm and the equal 
energy loss for both gluon and quark jets in central Au + Au collisions at -y/i = 200 
GeV. The numerical results are shown in Fig. ([6]) where the suppression factors are 
also included. The equal energy loss for both gluon and quark jets for (anti-)proton 
production in central Au + Au collisions results in the suppression factor R P AA rs 0.2, 
equal to R AA ~ 0.2 shown in Fig. ([5]). As shown in Fig. ([3]), high p T (anti-)proton 
production is dominated by gluon fragmentation. Compared with the case of parton 
energy loss AE g = \AE qi gluon jets with less energy loss AE g = AE q may have 
larger probability to pass through the hot medium and this gives a larger percentage of 
contribution for (anti-)proton production, so one can see in Fig. §6§ that (anti-)proton 
spectra with parton ernergy loss AE g = AE q are enhanced although the energy loss 
parameter e is adjusted large to 1.98 GeV/fm instead of 1.68 GeV/fm. 

Shown in Fig. ([7]) are p/p ratios in p + p and central Au + Au collisions at a/s = 
200 GeV. According to our NLO calculations on the fractions of contributions to (anti- 
proton from the fragmentaions of quark and gluon jets shown in Fig. ([3]), one know 
that the ratio of gluon to quark jet contributing to (anti-)proton production decreases 
with the (anti-)proton transverse momentum. Because most of antiprotons come from 
gluons while protons come from both valence quark and gluon fragmentation (31], the 
ratio of antiproton to proton production cross section should also decrease with their pt 
, as our NLO calculations show in Fig. (J7J). Similar to the LO calculations in Ref. [31J, 
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Figure 7. (Color online) The p/p ratios at mid-rapidity (\y\ < 0.5) as functions of pr 
in p + p (solid curve) and — 12% Au + Au collisions at y/s = 200 GeV. The data are 
from [321133]. 

the ratio will always be smaller than 1 because there are much more quark jets than 
anti-quark jets produced in hard scattering in p + p collisions due to baryon number 
conservation. In central Au + Au collisions the non-Abelian feature of parton energy 
loss [AE g = 9/AAE q ) reduces the contribution proportion from gluon jet fragmentations 
compared with that in p + p collisions, so the large transverse momentum p/p ratio in 
central Au + Au collisions is found to be smaller than that in p+p collisions. If gluon jets 
lose less energy [AE g = AE q instead of AE g = 9/4AE q ), the contribution from gluon 
jet fragmentations is enhanced, therefore there is a larger p/p ratio with AE g = AE q 
than that with AE g = 9/AAE q in central Au + Au collisions. In another point of view, 
equal energy loss for both quark and gluon jets causes to equal suppression for both 
p and p spectra (shown in Fig. (jEJ)), so the large transverse momentum p/p ratio with 
AE g = AE q in central Au + Au collisions is equal to that in p + p collisions. 

In Fig. (jHJ) we make a comparison of p{p)/tt + {tt~) ratios with 2 sets of parton energy 
loss (AEg = 9/AAE q and AE g = AE q ) in central Au + Au collsions at y/s = 200 GeV. 
On one hand, if gluon jets lose less energy (AE g = AE q instead of AE g = 9/4AE q ), 
more gluon jets are survived to contribute to hadron production when they encounter 
multiple scattering in the hot and dense matter created in central Au + Au collisions. 
The gluon contribution has an enhanced fraction while the quark contribution has a 
reduced fraction, and because high p T proton (anti-proton) production is dominated by 
gluon fragmentation whereas high px pion is dominated by quark fragmentation, the 
pip) i 'tt + (ti~) ratios with AE g = AE q are larger than those with AE g = 9/4AE q shown 
in Fig. (jHJ) . On the other hand, equal energy loss for both quark and gluon jets washes out 
the suppression difference between (anti-)proton and pion spectra, Raa(.Pt) ~ R\a{Pt)i 
shown in Fig. (jSJ) and ([6]). In fact, the p/rt ratio with parton energy loss AE g = AE q in 
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Figure 8. (Color online) The p/ir + and p/ir ratios at mid-rapidity(|y| < 0.5) as 
functions of pt in p + p collisions (solid line) and — 12% Au + Au collisions at y/s = 
200 GeV. The data are from [321E3]. 




Figure 9. Conversions between quark and gluon jets via both elastic gq(q) — + q(q)g 
and inelastic qq <-> gg scatterings with thermal quarks and gluons in the hot and dense 
matter. 



central Au + Au collisions can be written as 

/ P \AuAu i R P AA \fP\pv _ iP 



*7T R\a 71 71 



)(£1)pp « (JL)PP. (10) 



That's why our numerical results for large transverse momentum p/n ratios with 
AE g = AE q in central Au + Au collisions approach to those in p+p collisions, consistent 
with experiment findings[32, 33J shown in Fig. (jSJ). 



4. Scenario II: Strong Conversion between hard quark and gluon jets 

The conversion between quark and gluon jets is a mechanism introduced in Ref . [3U [35] to 
explain the discrepancy for p/w ratios between the experimental data and the theoretical 
calculations. Although a conversion enhancement factor K = 4 ~ 6 would be needed 
to explain experimental data observed by STAR and further theoretical studies will 
be needed to understand how we can get such a large conversion enhancement factor, 



A study on the anomaly of p over ir ratios in Au + Au collisions with jet quenching 13 





io- 


r 


I 


i i | 1 1 1 1 T= 


> 

ID 


10- 3 






STAR Au+Au 0-12% 
no -conversion £ =1.68 


o 


io- 4 






conversion £_=1.68 T 


N/dydp i 


10 5 
10 s 






f*f¥Tl V**I*ClflTl P — - 1 1 St — 


Q. 


10 7 




200GeV 




1/(271 


IO" 8 
10" 9 




Au+Au -> 0-12 




% n=1.2p T 

""Ts 




0.6 






• PHENIX 0-10% 1 


< 
< 
DC 


0.4 
0.2 



!i 








i i i i i i i i i i i i 






5 10 


15 20 



P T [GeV] 



Figure 10. (Color Online) The charged pion spectra and nuclear modification factors 
with and without conversion in central Au + Au collisions at y/s — 200 GeV. The data 
are from [331 Se- 
this mechanism indeed leads to an increase in the final number of gluon jets in central 
Au + Au collisions than the case without conversions [33j , and thus provide a promising 
way to explain the striking experimental data on proton over pion ratio observed at 
STAR collaboration [22]. In this section we will apply the jet conversion mechanism 
into the NLO pQCD parton model to study (anti-)proton and charged pion spectra and 
their ratios in heavy ion collisions. 

Conversions between quark and gluon jets could happen via both elastic gq(q) — > 
q(q)g and inelastic qq «-> gg scatterings with thermal quarks and gluons in the hot 
and dense matter as illustrated in Fig. (Q. For a parton jet created at a spatial point 
in Au + Au collisions region and escaping out off the dense matter in the transverse 
plane, multiple scattering with the medium particles wears off the jet energy and also 
stimulates the jet to convert to its partner. The conversion possibility for the jet passing 
through the matter is in turn connected with its energy and the medium particle density 
distributed along its passing path. 

As studied in Ref. [34j [50], the quark-to-gluon conversion possibility is found to 
be larger than the gluon-to-quark possibility at the same jet energy, and about 30% 
quarks for net quark-to-gluon jet conversions are needed to account for the observed 
ratios in the experiment in central Au + Au collisions. For simple evaluation for net 
quark-to-gluon conversions, we neglect gluon-to-quark jet conversions, and assume the 
quark-to-gluon jet conversion has an averaged possibility 45% all over its path. In our 
actual calculation, the quark and gluon jet contributions to the observed hadrons can 
be separated within a NLO pQCD parton model in A+A collisions. Once quark jets are 
created anywhere inside the medium, 45% quark jets are marked as gluon jets who will 
encounter multiple scattering and then fragment into hadrons. After a quark-to-gluon 
jet conversion happens at some point on the path in the medium, the newcomer (gluon) 
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Figure 11. (Color online) The separated p and p spectra and the suppression factors 
with and without conversion in central Au + Au collisions at y/s = 200 GeV. The data 
are from [33|47]. 




Figure 12. (Color online) The p/p ratio at mid-rapidity (\y\ < 0.5) as functions of pr 
in p + p (solid curve)and — 12% Au + Au (with and without conversion) collisions at 
^/s = 200 GeV. The data are from [32l [33] . 



will have energy loss AE g = 9/4AE q different from the energy loss AE q for its parent 
jet (quark). Because the "new" gluon is easier absorbed by the hot and dense matter 
than the parent quark, the total jet number will be reduced by net quark-to-gluon 
conversions, and there are actually about 30% of final surviving conversion jets. The 
ratio 30% is consistent to the studies in Ref. [34], [50j . In our calculations jet energy loss 
depends on the local parton density and total propagation length. A more sophisticated 
consideration will be needed in the further studies for the net quark-to-gluon conversion 
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Figure 13. (Color online) The p/tt + and p/k~ ratios at mid-rapidity < 0.5) as 
functions of pt in p + p and central Au + Au collisions with conversion (dotted-dashed 
line) and without conversion (dotted line) between quarks and gluons. Data points are 
taken from [32J, [33] . 

dependent on the local parton density and total propagation length. 

The appearance of net quark-to-gluon jet conversions will increase the total number 
of final gluon jets passing through the hot QGP medium. In such a conversion framework 
with gluon energy loss larger than quark energy loss, AE g = 9/AAE q , the highly hadron 
yield should be smaller than that without conversion. To fit data for hadron spectra, 
the energy loss parameter e should be adjusted to be smaller than 1.68 GeV/fm. 

Shown in Fig. ( flOi) are the charged pion spectra and nuclear modification factors 
with and without conversion in central Au + Au collisions at ^/s = 200 GeV. Our 
numerical results with conversion mechanism and eo = 1.18 GeV/fm are found to fit 
data well. With this energy loss parameter we then show in Fig. (fTTjl separately p and p 
spectra and the suppression factors with and without conversion in central Au + Au 
collisions. The spectra with conversion are found to fit data better than without 
conversion, and the nuclear modification factor for proton production with conversion in 
Fig. ( ITT]) is similar to that with conversion for pion production in Fig. ( fTU]) . R P AA ~ R\a 
due to the increasing number of converted gluon jets which dominate proton productions. 
However, Raa for anti-proton with conversion is slightly smaller than that for proton 
with conversion, Raa/^aa ~ 0.75, because most of antiprotons come from gluons with 
energy loss larger than quark while protons come from both valence quark and gluon 
fragmentation. 

In Fig. (Fl2|) we plot the p/p ratio as a function of Pt in p + p (solid curve) and 
central Au + Au (with and without conversion) collisions. Of interest is that the ratio 
p/p with conversion is found to be approximately equal to that without conversion, 
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and fit data better than p/p in the scenario of equal energy loss for both quark and 
gluon jets in central Au + Au collisions. In fact, according to the result in Fig. (fTTj) . 
R AA /R P AA ~ 0.75, one can obtain (j>/p)auAu ~ 0.75(p/p) pp , consistent with the numerical 
results in Fig. (fT2l) . The net quark-to-gluon jet conversions enhance the total number 
of the final gluon jets as compared with no conversion, and most of antiprotons come 
from gluons while protons come from both valence quark and gluon fragmentation, so 
the ratio p/p with conversion should be smaller than that with no conversion. But 
compared with the no conversion case with eo = 1.68 GeV/fm, the smaller energy loss 
parameter €q = 1.18 GeV/fm for the conversion case should give rise to the p/p curve 
for the conversion case a bit closer to the curve for the p + p case. Therefore, the ratio 
p/p with conversion is found to be approximately equal to that with no conversion in 
central Au + Au collisions. 

Shown in Fig. ( 1131) are the p/ix + and p/tx~ ratios in p + p and central Au + Au 
collisions with conversion (dotted-dashed line) and without conversion (dotted line) 
between quarks and gluons. The ratios of (anti)-proton over pion with conversion are 
found to fit data not bad. Net quark-to-gluon jet conversions effectively increase the 
total number of final gluon jets surviving in the fire ball. Therefore the final fractions 
of gluon and quark jets contributing to the hadron spectra will be modified by the net 
quark-to-gluon conversion. Furthermore, our numerical results in Fig. ([3]) show that 
even if without net quark-to-gluon conversions high pt proton (anti-proton) production 
is dominated by gluon fragmentations whereas high p^ pion is dominated by quark 
fragmentations in central Au + Au collisions. Consequently, larger p/7x ratios are found 
in central Au + Au collisions with conversion than without conversion, a°uAu > 
(^)Zir- and > (£)KT shown in Fig. (T3). Similar to the casein Fig.©, 

the ratios of (anti)-proton over pion in central Au + Au collisions with conversion are 
smaller than those in p + p collisions. According to R AA < R AA ~ R AA (= R-aa) from 
Fig. JE3KHH), one can get 

( A)^ u = ( |^ )( A r ^ ( A r (ii) 

1X^ R\a 71 71 

{ V_ )A uAu = (|M)(A r < (_L)» (12) 

TX R\ A TX TX 

as one can see in Fig. ( JT3l) . p/ix + ratios in central Au + Au collisions with conversion are 
similar to those in p + p collisions. However p/tx~ ratios in central Au + Au collisions 
with conversion are smaller than those in p + p collisions because the number of anti- 
quark jets is smaller than that of quark jets and the effect of conversion mechanism is 
not obvious for anti-quark jets. 

We note that in the conventional jet quenching theory, because a gluon jet always 
lose more energy than a quark jet with the same initial energy in a hot QCD medium 
(except in the case when both quark jet and gluon jet lose all of their energies), and 
at high pt vt meson is dominated by valence quark fragmentation whereas proton is 
dominated by gluon fragmentation, the p/7x ratio at large transverse momentum in 
p + p collision should be no less than those in A + A collisions. The two scenarios 
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considered in this paper will raise the p/ir ratios in A + A collisions therefore becoming 
closer to those in p + p, but will never go beyond. Only in the extreme case in A + A 
collisions at large Pt, all quark jets are converted into gluon jets, or equivalently there 
are only gluon jets, we can get higher p/7r ratio in A + A collisions at large pt- Though 
current data at STAR seems to imply that proton over pion ratio in Au + Au reactions 
is enhanced relative to that in p + p collisions, with large error bars at high p^, we still 
can not reach a unambiguous conclusion. If in the future more accurate experimental 
measurements confirm that at high p? region, p/it ration in A + A collision is enhanced 
significant compared with that in p+p collision, a dramatic change of our understanding 
of parton energy loss mechanism will be needed. Because we focus on p/n ratios at 
the high px region, in our calculations we have ignored the contribution of parton 
recombination/coalescence mechanism in the QGP. We expect this negligence underlines 
the rather large discrepancy between our numerical simulations and experimental data 
in the intermediate px region, where the recombination/coalescence mechanism in the 
hot QGP should play a very important role [SH E2J [53], [5U 155] . 

5. Summary and discussions 

In this paper, we investigate p/n + and p/ir~ ratios at large transverse momentum in 
Au + Au collisions with jet quenching based on an next-to-leading order pQCD parton 
model. Two scenarios are considered to explain the anomaly of p/n ratios observed 
by STAR Collaboration that p/rr ratios at high p? in central Au + Au collisions [32] 
approach those in p+p and d+Au collisions [33] . Firstly, we investigate the consequence 
in Scenario I where gluons lose the same amount of energy in QGP as quarks and 
demonstrate that it will enhance thep/7r + andp/7r~ ratios for central Au + Au collisions. 
Secondly, we explore the effect of strong jet conversion in Scenario II where a net of 
quark-to-gluon conversion also result in similar p/n + andp/7r - ratios in central Au + Au 
and p+p collisions. In both scenarios, the final number of gluon jets passing through the 
hot and dense medium is increased as compared to the original jet quenching picture, 
and jet conversion scenario is favorable when p/p ratio is considered because p/p in the 
latter scenario is found to fit data better than that in the former scenario. This kind of 
studies should improve our understanding on the flavor dependence of energy loss and 
lead us to a unified picture of jet quenching. 

We note that in our approach the effect of energy loss is taken into account through 
the medium-modified effective fragmentation functions and the evolution of the hot 
nuclear medium is described by a simple 1+1D Bjorken evolution based on Hard- 
Sphere geometry. This approach has given very well descriptions on inclusive pion 
suppression jS], di-hadron production [37] and photon-tagged hadron production |56j, 
but several improvements will be needed to make a more realistic consideration of parton 
energy loss in the hot medium created at RHIC. One direction is to simulate the medium 
evolution with the 3+1D hydrodynamics model as in Ref. [571 EH], which may give 
corrections to the total parton energy loss though the effect of local density dropping due 
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to the transverse expansion of the medium will be offset by the longer path length in the 
transverse direction [59] in the QGP phasdj]. Also one could study the in-medium parton 
shower [57J[6lJ[62] to make a more differential study of jet quenching than the adoption 
of the effectively modified parton fragmentation function. Furthermore a more detailed 
medium dependence of jet energy loss can be obtained by introducing the energy loss 
probability distributions [53, [581 E31 Ell E2], which will reduce the difference of energy 
losses of quark and gluon jets by properly incorporating kinematical constraints, and 
thus raise the p/n ratios. Last but not the least, in our calculation we have assumed 
that the pions and protons are formed outside of the medium, and the hadronization of 
pions and protons in A+A collisions are the same as that in p+p collisions. Larger p/it 
ratios may be obtained and help to resolve the p/it puzzles observed in A+A collisions at 
STAR if one assumes that protons in A+A collisions can be formed inside the medium 
due to its larger mass relative to pion , though so far how to model hadronizations of 
partons with large pt in medium and whether it is distinct from that in vacuum is still 
an open question and under hot debate. Nevertheless, intensive theoretical studies will 
be needed to see whether we can resolve the puzzle of p/n ratios in a more beautiful 
and natural way by improving our model and including all the relevant effects discussed 
above. 
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